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ABSTRACT 



We discuss the consequences of almost rectilinear acceleration of protons to ex- 
tremely high energies in a reconnection region on the surface of an accretion disk which 
surrounds a central black hole in an active galaxy. The protons produce 7-rays and 
neutrinos in interactions with the disk radiation as considered in several previous pa- 
pers. However, in this model the secondary 7-rays can initiate cascades in the magnetic 
and radiation fields above the disk. We compute the spectra of 7-rays and neutrinos 
emerging from regions close to the disk surface. Depending on the parameters of the 
reconnection regions, this model predicts the appearance of 7-ray and neutrino flares 
if protons takes most of the energy from the reconnection region. In contrast, if leptons 
take most of the energy, they produce pure 7-ray flares. The 7-ray spectrum expected 
in the case of hadronic cascading is compared with the spectrum observed during the 
flare in June 1991 from 3C 279. The neutrino flares which should accompany these 
gamma-ray flares may be detected by future large scale neutrino telescopes sensitive 
at - lO'^ TeV. 
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1 INTRODUCTION 

Flares on very short time scales in TeV 7-rays have been re- 
cently detected from Mrk 421, Mrk 501, and 1ES2344-I-514 
(e.g. Buckley et al. 1996, Gaidos et al. 1996, Aharonian et 
al. 1997, Catanese et al. 1998) and in the EGRET en- 
ergy range from many other blazars (e.g. von Montigny et 
al. 1995, Mattox et al. 1997). Flares observed with a ~ 1 day 
time scale are usually interpreted in terms of a relativis- 
tic shock moving through the jet with a Lorentz factor of 
the order of ~ 10 (e.g. Maraschi, Ghisellini & Celloti 1992, 
Mannheim & Biermann 1992, Dermer & Schlickeiser 1993, 
Bednarek 1993, Sikora, Begelman & Rees 1994, Blandford & 
Levinson 1995, Ghisellini & Madau 1996, Protheroe 1997). 
However, the production of such flares may also be a natu- 
ral consequence of energization of particles in a way sim- 
ilar to that causing solar flares, but in this case occur- 
ring in regions of magnetic reconnection which may exist 
on the disk surface (Haswell, Tajima & Sakai 1992, Lesch 
& Pohl 1992). Reconnection of magnetic field may also oc- 
cur in the jet (e.g. Romanova & Lovelace 1992, Bisnovatyi- 
Kogan & Lovelace 1995). The processes of 7-ray produc- 
tion by leptons accelerated in electric fields generated in the 
jet has been recently discussed by Bednarek & Kirk (1995, 



henceforth BK95) and Bednarek, Kirk & Mastichiadis (1996, 
henceforth BKM96). 

The acceleration of particles by almost rectilinear elec- 
tric fields has some advantages in comparison to the stochas- 
tic acceleration by a relativistic shock since, in principle, 
particles can reach very high energies on shorter time scales. 
However, the process of reconnection of magnetic fields in 
astrophysical environments is poorly understood. Recently, 
Haswell, Tajima & Sakai (1992) have developed a model 
of explosive reconnection of magnetic fields on the sur- 
face of an accretion disk. They argued that particles (elec- 
trons, protons) can be accelerated to extremely high ener- 
gies (~ ]^q2o-2i qY). The consequences of this model for 
7-ray production in blazars has been investigated recently 
by Bednarek (1997, henceforth B97), assuming that acceler- 
ated leptons develop cascades in the magnetic field through 
production of 7-rays by the synchrotron process, and by the 
creation of secondary pairs by 7-rays in the magnetic field. 
The 7-ray spectra, emerging from a region on the inner part 
of an accretion disk around a massive black hole, are in 
agreement with observations of blazars. 

In this paper, we discuss the consequences of accelera- 
tion of protons in almost uniform electric fields induced by 
magnetic reconnection on the disk surface. Protons acceler- 
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ated to extremely high energies produce 7-rays and neutri- 
nos through the decay of pions produced in collisions with 
disk photons. Neutrinos freely escape, but 7-ray photons 
initiate cascades in the magnetic and radiation fields above 
the disk (B97). The conditions for acceleration of protons, 
and their propagation in the disk radiation, are discussed in 
Sect. 2.1 and 2.2. We compute the spectra of neutrinos and 
7-rays (Sect. 2.3), and describe the details of the cascades 
initiated by these 7-rays (Sect. 3). Example 7-ray spectra 
are compared with observations of the June 1991 flare from 
3C 279, and a prediction of the neutrino emission accompa- 
nying this flare is given in Sect. 4. 



2 ACCELERATION OF PROTONS DURING 
EXPLOSIVE RECONNECTION IN THE 
INNER PART OF AN ACCRETION DISK 

We make use of the model of Haswell et al. (1992) in which 
particles can be accelerated to extremely high energies dur- 
ing explosive reconnection of the magnetic field on the sur- 
face of an accretion disk in the central parts of active galac- 
tic nuclei (AGN). Let us assume that reconnection of the 
magnetic field occurs over a distance /rec somewhere in the 
inner part of the optically thick accretion disk. The local 
radiation field in this region can be described by a black 
body spectrum with temperature T. The induced electric 
field strength Ef can be estimated as Ei = cB, where B 
is the reconnecting magnetic field strength and c is the ve- 
locity of light. Then, the maximum possible energies, which 
can be reached by particles (in eV), are 



; 3 X W^^Bshi, 



(1) 



where B = 10^ B5 G, 



10"Zii cm. 



2.1 Conditions for acceleration of protons 

Protons accelerated in the electric field lose energy by 
pair and pion production in collisions with the black body 
photons after reaching the thresholds for these reactions. 
Since the threshold for pair production is about two or- 
ders of magnitude lower than that for pion production, 
pairs are first injected into the electric field in the reconnec- 
tion region. The characteristic distance scale Ap.^g± (in cm), 
for this process to occur efficiently, can be estimated by 

\7e± ~ (o-p^^e±'^bb)"^ ~ 5.1 X (2) 

where ap^^^± « 10"^'' cm'^ is the cross section for 
pair production in collisions of protons with photons (Max- 
imon 1968, Chodorowski, Zdziarski & Sikora 1992), ribb is 
density of the black body photons, and the local disk tem- 
perature in the reconnection region is expressed in units of 
T4 — T/IO* K. We shall show later that this process is not 
an important energy loss mechanism for accelerated protons 
since its inelasticity coefficient is small (~ 2 x 10~^). Never- 
theless, the injected secondary pairs can initiate inverse 
Compton pair cascades in the disk radiation and the elec- 
tric field of the reconnection region (BK95). The products 
of such a cascade can saturate the external electric field if 
the length of the reconnection region is longer than (BK95) 

A,at ~ lOOAtrip ~ 3.85 x 10"T4~^ (3) 



where Atrip ~ (o"trip"-bb)~^ is characteristic distance scale for 
the triplet pair production by electrons in collision with 
soft photons, and atrip ~ (tt/50 (Mastichiadis 1991), (tt is 
the Thomson cross section. In this case, most of the energy 
extracted from the electric field is transferred to secondary 
pairs if 

^rcc > Asat. (4) 

Eq. ^ can be written as a condition on the temperature of 
the black body radiation in the reconnection region such 
that the electric field is saturated, 

T4 > r4,tr ~ l5.7^~'''^ (5) 

If the temperature of the black body radiation is lower 
than r4,tr then the reconnection energy can efficiently accel- 
erate relativistic protons. These protons can then produce 
pions in interactions with the disk radiation both during 
the acceleration stage and after emerging from the reconnec- 
tion region, provided that the electric field is strong enough 
to allow them to reach energies, ii'p,min (in eV), above the 
threshold for pion production. 



Ep > Ep 

,min 



(6) 



The characteristic distance scale for pion production by pro- 
tons can be estimated by, 

ApT^ ~ (o-p^^ribfc)"^ ^ 1.7 X 10^''r4~^ (7) 

where apj,r ~ 3 x 10~^* cm~'^ is the cross section for 
pion production in collisions of protons with photons 
(Stecker 1968). Comparison of Ap-Y,r with the characteristic 
distance scale of the reconnection region, which is the inner 
radius of an accretion disk r-m, allows us to estimate the pa- 
rameters of the accretion disk for which the proton energy 
losses by pion production become significant. We find that 
the minimum local disk temperature in the reconnection re- 
gion depends on the characteristic disk dimension, 

T4 > 1.2r;^'^\ (8) 

where nn = 10^'*ri4 cm. 

The maximum proton energies, estimated using Eq. 1, 
can be expressed in terms of the temperature of the black 
body radiation in the reconnection region if we use the com- 
mon assumption of equipartition of the magnetic and radi- 
ation energy densities close to the accretion disk. However, 
we note that Haswell, Tajima & Sakai (1992) argue that the 
magnetic field in the inner disk can be significantly amplified 
as a result of its differential rotation from the values close 
to the equipartition (Bcq ~ 4OT4 ) up to values of the order 
of -B ~ 300-Bcq. In such a case, the formula for the proton 



energy (Eq. 1) can be rewritten in the form 
Ep ^ 3.6 X lO^Tihi. 



(9) 



Since efficient proton acceleration requires the electric field 
not to be saturated by pairs, the disk temperature in the 
reconnection region cannot exceed that given by Eq. ^ and 
we can express the maximum possible proton energy (in eV) 
in terms of the local disk temperature 



1.4 X 10'^^ Tr^. 



(10) 



The assumption of equipartition between the primary mag- 
netic field energy density and the radiation energy density 
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Figure 1. The conditions for cascading initiated in the radiation 
field of the accretion disk by hadrons accelerated on the disk sur- 
face. For the parameters in the areas marked by 'hadrons escape', 
the protons escape from the radiation of the disk without signifi- 
cant energy losses. In the area 'no e+e^' they have energies below 
the threshold for pair production and in the area 'e+e~', they 
inject pairs. For parameters in the area marked by 'hadronic 
cascades', protons efficiently lose energy on pion production and 
the pion decay products initiate cascades in the magnetic and 
radiation fields. In the area 'leptonic cascades', the energy of the 
electric field is extracted mainly by secondary leptons which ini- 
tiate cascades in the magnetic and radiation fields (Eq. 5). The 
thick long-dashed line shows the limit for efficient extraction of 
energy from hadrons (Eq. 8) for a disk inner radius of Tin = 10^^ 
cm, and the thin long-dashed lines give the limits for rin = 10^* 
cm (right) and 10^^ cm (left). 



i 




Figure 2. Schematic representation of the scenario discussed in 
this paper. Relativistic protons with Lorentz factors 7p are in- 
jected at a reconnection region which is at a distance ri^j from 
the center of the disk. A proton interacts at height Xp with a soft 
photon from the disk, having energy e, emitted at point A on the 
disk defined by the angles a and 0. The energy t is determined by 
the local disk temperature T(r) at point A. As a result, pions are 
produced which decay into 7-rays, pairs, and neutrinos. The 
surviving proton, with the Lorentz factor 7p/ , may again interact 
at some distance from the disk. 



allows us the derive the minimum length (in cm) of the re- 
connection region for which accelerated protons can reach 
energies above the thresholds for pair production and 
pion production (Eq. |) giving 

Itt > 6.9 X 10■^^4"^ (11) 
and 

Z^cc > 8.3 X 10''^4"^ (12) 

in cm, respectively. The conditions for acceleration and 
propagation of protons (Eqs. 5, 8, 11, and 12) allow us to 
distinguish three interesting scenarios which can be sepa- 
rated in the parameter space (shown in Fig. |l|) defined by 
the length of the acceleration region Ir^c and the temper- 
ature of the local black body radiation T. For the regions 
labeled by 'hadrons escape', the proton energy losses during 
acceleration and subsequent propagation in the disk radia- 
tion are negligible, allowing protons to leave the accretion 
disk with very high energies. If the acceleration distance is 
high enough, the protons can reach the threshold for pair 
production, injecting a few secondary pairs (area marked by 
'e^e~', between the dot-dashed line (Eq. [13), dotted line 
(Eq. ^ and long-dashed line (Eq. ^). Below the dotted 
line, protons do not inject pairs (area marked by 'no e"'"e~ '). 
If protons reach the threshold for pion production (Eq. h2l 
marked by the dot-dashed line), they suffer strong energy 
losses. The proton energy losses becomes significant for disk 



temperatures above the values estimated by Eq. |8| (marked 
by the long-dashed lines for different disk inner radii). In 
the area marked by 'hadronic cascades', the initial proton's 
energy is transferred to pions whose decay products (lep- 
tons, photons) initiate cascades in the magnetic and radi- 
ation fields above the disk. In the present paper, we shall 
concentrate on this case. In the last area, marked by 'lep- 
tonic cascades', the production of pairs by protons and 
subsequently by these pairs in the inverse Compton pair 
cascade developing in the electric field, is very efficient, even 
during the acceleration stage in the reconnection region. 
These secondary pairs may even saturate the electric 
field (BK95). Most of the energy of the electric field is then 
transferred to leptons (and 7-rays) which initiate cascades 
in the magnetic and radiation fields above the disk identical 
to those initiated by primary leptons (B97). Note that the 
processes occurring for the parameters in the area 'hadronic 
cascades' turn out to be the production of 7-ray and very 
high energy neutrino fiares, whereas for the processes oc- 
curring for parameters in the area 'leptonic cascades' only 
7-ray flares are produced. 

2.2 Propagation of relativistic protons in the disk 
radiation 

The protons, escaping from the reconnection region with 
energies estimated by Eq. (ffl or (llOh, produce pions in col- 
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lisions with the disk photons on a characteristic distance 
scale given by Eq. (Q). We simulate propagation of protons 
through the anisotropic disk radiation by using the Monte 
Carlo method. It is assumed that protons are injected on the 
surface of a thin disk at the distance rinj from its center (see 
Fig. We consider a simple disk model in which the disk 
emits black body radiation with characteristic temperature 
depending on the distance from the disk center according 
to the model of Shakura & Sunyaev (1973). Therefore, the 
disk radiation is defined by the disk inner radius rin and its 
temperature Tin at rin. 

In order to simulate the propagation of protons and the 
production of 7-rays and neutrinos, we first determine the 
point of interaction of the proton with a soft photon coming 
from the disk. The interaction length, Ap-y-n-, for the proton 
in the disk radiation, as a function of its energy Ep (Lorentz 
factor 7p and velocity /3p) and its distance, x, from the disk 
surface is calculated, by integrating the formula 



n(e, n)ap-fTr{e') — dedfl, 



(13) 



where n(e, f2) is the differential density of photons emit- 
ted by the disk, (Jp^.,r{e') is the energy dependent cross 
section for pion production in proton-photon collision 
(Stecker 1968), dQ = dcjidfi is the solid angle of incident 
disk photons, /x = cosa (see Fig. e = e~fp{l + f3p^) and e 
are energies of soft photons in the proton and LAB frames, 
respectively. To find out if protons have a chance to interact 
with the disk photons for specific disk parameters, the op- 
tical depth for protons is obtained by integrating along the 
proton path 



\^AEp,x)dx. 



(14) 



The optical depth for protons is shown in Fig. |3|a as a func- 
tion of proton energy, and for different parameters of the disk 
radiation, rin and Tin. It is clear that for the disk parameters 
expected in AGNs protons have a significant probability of 
interacting, and may even undergo multiple interactions. In 
order to find out if the competing process of pair pro- 
duction may become an important energy loss mechanism 
for protons, we calculate the energy loss rates in the disk 
radiation for pair production at a fixed distance, x, from 
the disk surface 



dEpjx) 
dx 



n{e, D,)ap^^± {e)K{e') — dedD,, 



(15) 



where the cross section for pair production, o-p^^±{e), 
and the inelasticity K{e) are taken from Chodorowski, 
Zdziarski & Sikora (1992). The optical depth for the total 
energy loss by protons during their propagation through the 
disk radiation from the disk surface to infinity is 



^-^Ux. 

dx Ep 



(16) 



This value is plotted in Fig. ^b, for the same disk param- 
eters as used in Fig. ^a. It is evident that proton energy 
losses on pair production are not important except for 
the case of very large disks with high inner disk tempera- 
tures (rin > 10^^ cm and Tn > 10^ K). Such disks, with 
luminosities Tdisk > 7 x lO^'' erg s~^, have not been ob- 
served. The blazar 3C 273 has a disk luminosity a factor of 



about 3 lower, Tdisk ~ 2 x lO"** erg s~^ (Lichti et al. 1995), 
but this seems to be exceptionally high. Therefore, we ne- 
glect the infiuence of pair production on the propagation 
of protons since these losses can be neglected in comparison 
to pion production losses at energies above i?p,min (Eq. 
However this process can become important if proton ener- 
gies drop below the threshold for pion production far from 
the accretion disk. Such pairs could contribute to the forma- 
tion of a broad synchrotron bump at low energies observed 
in blazars. 

If protons are very close to the disk, their interaction 
length does not change significantly with distance xo from 
the disk. Using the Monte Carlo method, the point of inter- 
action of a proton with the disk photon can be then deter- 
mined from 



Xp = xo " Xpj.^ln{l ~ Pi), 



(17) 



where Pi is a random number. However at larger distances 
from the disk the radiation becomes anisotropic, and the 
proton interaction length can change significantly. Then the 
distance of proton interaction, Xp, is obtained by solving 



Pi = exp[- 



Ap-)/7r(/) dl\ 



(18) 



where I is the distance measured along the proton's path. 
However, if 



Pi < exp[- 



Ap"j'7r(/) dl] 



(19) 



the proton is assumed to have escaped from the disk radia- 
tion. The condition applied in order to distinguish between 
quasi-isotropic and anisotropic cases is this same as used in 
Bednarek (B97) (see footnote on page 5 in that paper). 

Having selected randomly the distance, Xp, at which 
proton interacts with a disk photon, we have to sample the 
cosine angle, fi, of the photon from the disk which interacts 
with the proton. The possible range of photon directions 
(from /i = — 1 to /imax, which corresponds to the outer disk 
radius) is divided into one hundred intervals with width A/x, 
corresponding to annular rings on the accretion disk cen- 
tered on Hnj. The ring defined by fj,ph{Ep,Xp), from which a 
disk photon is interacting with the proton, is selected ran- 
domly by solving 



P2 = 



Mph(-Bp,a:p) ^ 

J2 f f n{e,D,)ap.,Tr{(^')^ded(j)A^i 
-1 

Mmax ^ 

J2 1 1 n{e,Q,)ap^^{e')^ded(l)An 



(20) 



where P2 is a random number. 

Finally we have to sample the energy of the disk photon 
coming from the ring /Xph for a proton located at a distance 
Xp. To save computing time we approximate the black body 
spectrum of photons emitted at distance r from the disk 
center by monoenergetic photons with energies correspond- 
ing to the mean photon energy of the black body spectrum, 
e = 2.7kBT{r), where fee is the Boltzman constant, and 
T(r) = Tin(r/rin)~"^^*, as predicted by Shakura & Sunyaev 
model. In such a case, the energy of photons coming from 
a specific ring Uph is uniquely defined by the azimuthal an- 
gle cj) (see Fig. bh. We can sample the angle cj) randomly by 
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Figure 3. (a) Optical depth for pion production collisions of protons with the disk photons; (b) optical depth for energy loss of protons 
by production of pairs during propagation. The full curves show the results for disk inner radius rj^ = 10^'"' cm and different inner 
disk temperatures: Tin = 10* (bottom), 3 X 10*, and 10^ K (top); the dashed curves are for Tin = 3 X 10* K and different inner disk 
radii: rin = 10^* (bottom), and 10^® cm (top). 



solving 



Ps = 



TT 



(21) 



where P;i is a random number and (pmin is equal to 0, or 
is calculated for the given inner disk radius rin if the ring 
defined by n lies partially inside rin- The interacting pho- 
ton is now completely specified, allowing us to compute its 
energy in the proton frame e' = e(0e)7p(l + PpHph), model 
the interaction, and obtain the energies of pions produced 
in the proton frame. 

The above procedure is repeated until the relativistic 
proton escapes from the disk radiation. The parameters of 
produced pions are stored for further computation of neu- 
trino and 7-ray spectra. 

2.3 Spectra of neutrinos and gamma-rays from 
decay of pions 

Because of the complexity of the present model, in particular 

the cascading in the radiation and magnetic fields, we have 
had to make a number of approximations in order to save 
computer time. In particular, in the calculation of spectra 
of secondary particles produced in interactions of protons 
with soft photons was, of necessity, subject to approxima- 
tion. However, we include single pion production (in the A 
resonance region) although in an approximate way, and the 
effects of multiple pion production far away from the thresh- 
old is treated using an average multiplicity relation similar 
to that found in low energy pp interactions, 

< >oc (22) 

where s = rUp + 2eEp{l + (ipjj) is the center of momentum 
frame energy squared (in this section we define c= 1), and 
rup is the proton rest mass. The ratio of charged pions to all 



pions produced in a single multiple pion interaction is taken 

to be n^± jn-r: « 2/3. 

In the case of single pion production (A resonance de- 
cay) the pion energies in the center of momentum system 
(CMS) are given by, 

j^CMS _ 2mpe' -I- mi 



2^ 



(23) 



It is assumed that the A resonance decays isotropically, so 
that pions are produced isotropically in the CMS and we 
may sample their direction cosine angle using 



2P4 



1, 



(24) 



where P4 is a random number. The energy of the pion in the 
LAB frame is then 

Ei^"" = 7CMs(£;?^' + /JcmsP^'mx), (25) 

where 7cms = {Ep + e)/^/s, /3cms = ^1 - 1/7cmS' 
p^^^ is the pion momentum in the CMS. 

The pions decay isotropically in their rest frame pro- 
ducing: 7-rays, muon neutrinos, and muons with energies 
which we approximate by 

EJI = 0.5m,r, 



= 0.5{mi - m'l)lmn, 

and 

EZ = 0.5(m^ -I- mt,)/m^, 



(26) 
(27) 



(28) 

where mn is the rest meiss of a pion (we neglect the mass 
difference between neutral and charged pions). The cosine 
angles of decay products of pions can be sampled randomly 
from 

111, = 2P5 - 1, (29) 

where P5 is a random number, and Hi^i = —fiv. Their energies 

in the LAB frame are given by 



(30) 
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Figure 4. The spectra of 7-rays (full curves) and muon-neutrinos (long-dashed curves) from decay of pions which are produced in 
collisions of protons with the disk photons. The parameters used are: (a) Tin = 3 X 10" K, Ircc = 10^2 (giving Ep = 3 X lO^^ MeV); 



(b) Ti, 



10* K, Ir 



- 3 X 10 cm (Ep = 10 MeV). In each case results are given for r'in = 10 cm (bottom), 10 cm, and 10 cm 



(top). The dotted and dot-dashed curves in (a) show the spectra of neutrinos and 7-rays for protons with Ep = 3 X 10 MeV. 



where 7.^ — E^^^ /rriTr, and P-^ = y^l — 1/7^. The energies 
of muon-neutrinos from the decay of secondary muons are 
simply estimated by taking E[!; ~ 

The spectra of muon-neutrinos and 7-rays, averaged 
over the propagation of protons along the direction perpen- 
dicular to the disk surface, are shown in Fig. 4. Fig. ^ shows 
the spectra for the local disk temperature in the region of 
magnetic reconnection T = 3 x lO^K, /rcc = 10^^ cm, and 
T'inj = 1.5rin. For these parameters, protons are injected from 
the reconnection region with energies Ep = 5 x 10'^'^ MeV 
(see Eq. The full curves show the spectra of 7-rays for 
the inner disk radii: nn ~ 10^*, 10^^, and 10^^ cm. The corre- 
sponding neutrino spectra are shown by the dashed curves. 
The spectra of 7-rays and neutrinos in the case of protons 
injected with energies = 3 X 10^° MeV are shown by the 
dot-dashed and dotted curves, respectively. In Fig. Mb we 
show the spectra of neutrinos and 7-rays for T = 10* K, 
cm, Tinj — LS^in- For these parameters the 
proton initial energy is E^ — 10^^ MeV. 



3 CASCADE INITIATED BY SECONDARY 
GAMMA-RAYS 

As shown in Figs. 4, the 7-rays originating from decay of 
neutral pions have extremely high energies. Depending on 
the conditions above the reconnection region, they can de- 
velop cascades in the perpendicular component of the mag- 
netic field via magnetic pair production, and subsequent 
emission of synchrotron photons, provided that their ener- 
gies fulfill the condition (Erber 1966, B97) 



_(m,cC Be 



> 0.05, 



(31) 



where B± is the component of the magnetic field which is 
perpendicular to the photon propagation, E-y is the photon 
energy, Bcr = mlc^{eh)'^ = 4.414 x 10^^ G is the critical 
magnetic field strength, and rric is the electron rest mass. 
For hot disks having weak magnetic fields, the pair 



energy losses by ICS dominate over synchrotron losses, and 
the secondary pairs initiate ICS pair cascades (e.g. 
Protheroe, Mastichiadis & Dermer 1992, Coppi, Kartje & 
Konigl 1993). In the present work we discuss the case 
where the magnetic and disk radiation fields are initially 
in equipartition close to the surface. Therefore the mag- 
netic energy density dominates over the radiation energy 
density above the disk since the magnetic field drops more 
slowly with the distance from the disk than the radiation 
field (Bednarek, Kirk & Mastichiadis 1996b). 

The energies of 7-ray photons produced by protons and 
their distances from the disk on production are stored and 
used as input to the subsequent cascading in the magnetic 
field calculated by using the numerical cascade code de- 
scribed in (B97) . A sequence of different processes is involved 
during propagation of these 7-rays. First, the 7-rays cascade 
in the perpendicular component of the magnetic field above 
the disk, and we compute the spectra of secondary 7-rays 
and pairs emerging from this synchrotron cascade. 
At distances further from the disk, where the cascade in the 
magnetic field becomes inefficient, the propagation of 7-rays 
in the disk radiation is considered. If the inner disk temper- 
ature is high enough, some of secondary 7-rays are absorbed 
creating a second generation of pairs. These pairs cool 
predominantly in the local magnetic field by the synchrotron 
process, emitting a second generation of 7-rays. Some of the 
most energetic 7-rays from second generation can further be 
absorbed in the disk radiation producing tertiary pairs 
which cool by the synchrotron process as well. Since the opti- 
cal depth for tertiary 7-rays in the disk radiation is low, they 
usually escape unabsorbed. The emerging 7-ray spectrum is 
obtained by summing the contributions from all these pro- 
cesses described above. In the present work we use a more 
general prescription for the dependence of the perpendicu- 
lar component of the magnetic field along the jet {B±) than 
applied by Bednarek (B97). Here we use 



B±{x) = Boix/xiY 



(32) 
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Figure 5. The 7-ray spectra obtained in a sequence of processes, discussed in tlie text, wliich are initiated by primary 7-rays produced 
for the parameters of the model as in Figs. 4(a) and 4(b), respectively, and = 10^^ cm. The 7-ray spectra from decay of neutral pions 
produced by protons in collisions with the disk photons are shown by the dotted curves; the 7-ray spectra emerging from the cascade in 
the magnetic field above the disk are shown by the dashed curves; and the 7-ray spectra from absorption of the cascade 7-ray spectra in 
the disk radiation, and subsequent synchrotron cooling of secondary pairs are shown by the full curves. 



which for the cases with (5 < 1.5 guarantees that the 
magnetic energy density always dominates over the disk, 
radiation energy density above the disk provided that it 
dominates at the base of the jet (Bednarek, Kirk & Mas- 
tichiadis 1996b). 

The spectra of 7-rays from the synchrotron cascade, 
and the final 7-ray spectra escaping from the radiation field 
of the accretion disk, are shown in Fig. 5 by the dashed and 
full curves, respectively. They correspond to the initial 7- 
ray spectra from the decay of pions shown in Fig. 4, and for 
a disk inner radius rin = 10^^ cm. In these simulations, we 
apply the following parameters which define the structure of 
the magnetic field in the jet above the reconnection region: 
Bq is taken to be equal to the value of the magnetic field close 
to the disk surface, x\ = nnj, and /3 = 1.25. Our simulations 
show that the shape of the final 7-ray spectrum is not very 
sensitive to (3 for values in the range between 1. and 1.5. 



4 JUNE 1991 GAMMA-RAY FLARE FROM 3C 
279 AND EXPECTED NEUTRINOS 

The 7-ray spectra produced in the sequence of processes dis- 
cussed above are consistent with the spectra observed from 
FSRQ blazars. As an example, in Fig. ^ we compare the 
7-ray spectrum, shown in Fig. Bp,, with the observations of 
the flare from 3C 279 during June 1991 by the COMPTEL, 
EGRET and Ginga detectors (Williams et al. 1995, Knif- 
fen et al. 1993, Hartman et al. 1996). The accompanying 
neutrino flare (marked by the dashed curve) corresponds to 
~ 28 muon neutrinos per day per km^ with a typical en- 
ergy of ~ 10* GeV. Such a flux, although small, may still be 
detectable since it is many orders of magnitudes above the 
neutrino background induced in the atmosphere within the 
1° (Lipari 1993) shown by the dot-dashed curve in Fig. ^ 
The detection probability of neutrino telescopes operating 
at this energy can be estimated to be ~ 3 x 10"'^ (Gaisser, 



Halzen & Stanev 1995), giving a rate of ~ 0.1 muons per 
day per km^. Therefore, for typical flares occurring on a 
time scale of a week one coincident neutrino event should 
be detected provided that the neutrino is massless. Dur- 
ing a period of years, a large neutrino detector should be 
able to deflnitively answer the question of the importance 
of hadronic processes by measurements of correlations be- 
tween neutrino arrival directions and directions to known 
FSRQ blazars. 

The spectrum of neutrinos expected in this model is lim- 
ited to a relatively narrow energy range if compared to other 
hadronic models (e.g. Protheroe 1997, Nellen, Mannheim & 
Biermann 1993, Mannheim 1995, and comparison of differ- 
ent models in Protheroe 1998). The reason is that we have 
considered the case in which all hadrons are accelerated in 
the reconnection region to this same maximum energy, in 
contrast to other models mentioned in the Section 1, which 
assumed the injection of a power-law spectrum of protons. 
Therefore, the shape of the neutrino spectrum may also give 
us important information on the acceleration mechanism of 
hadrons in blazars. 



5 DISCUSSION AND CONCLUSION 

Some models of magnetized accretion disks predict the ap- 
pearance of relatively small, but efficient, acceleration re- 
gions close to the surface of the inner part of the accretion 
disk as a result of magnetic reconnection (Haswell, Tajima 

6 Sakai 1992, Lesch & Pohl 1992). In this paper, we discuss 
the consequences of acceleration of hadrons in such regions. 
Depending on the conditions in the acceleration region, and 
above the disk (see Fig. |l|), accelerated hadrons may either 
escape, initiate cascades after emerging from the accelera- 
tion region, or produce copious secondary leptons during the 
acceleration process which initiate subsequent leptonic cas- 
cading. Primary leptons can also initiate cascades, and this 
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Figure 6. The 7-ray spectrum obtained from the present model 
(full curve) is compared with observations of the June 1991 flare 
from 3C 279 detected by COMPTEL (C) (Williams et al. 1995), 
EGRET (E) (Kniffen et al. 1993) and Ginga (G) (Hartman 
et al. 1996). The parameters of the fit are as in Fig. The 
dashed curve shows the expected muon neutrino flux during this 
flare. The dot-dashed curves show the atmospheric neutrino back- 
ground (horizontal - upper curve, and vertical - lower curve) for 
a neutrino detector with 1° angular resolution . 



possibility has been already discussed (B97). In the analysis 
of hadronic cascades we included the processes of pion pro- 
duction by hadrons, synchrotron radiation of secondary 
pairs in the quantum domain, magnetic pair production 
by 7-ray photons and absorption of the 7-rays in the disk ra- 
diation. The 7-ray spectra emerging from such a sequence of 
processes resemble the spectra observed from FSRQ blazars 
(Figs 5). Their comparison with the June 1991 flare from 3C 
279 allows one to predict the strength of the accompanying 
neutrino flare (see Fig. ^). 

In the present model we have neglected the contribution 
of interactions of accelerated hadrons with the background 
matter, either during the acceleration stage, or after propa- 
gation through the disk radiation. Interactions of relativistic 
protons with the disk radiation dominate over the interac- 
tions with matter if the density of matter riH is lower than 



riH 



2 X 10"T4 particles cm"^. 



(33) 



If the density of matter is relatively high at larger distances 
from the disk, then protons, at this stage having lost most 
of their energy during propagation through the disk radia- 
tion, can also produce neutrinos by interactions with mat- 
ter. These neutrinos are typically of lower energies, and their 
power will be much less than in the higher energy neutrinos 
from interactions with radiation. 

Here we have considered only a single acceleration re- 
gion in order not to complicate the picture too much. How- 
ever, in principle many acceleration regions with different 
parameters could be present simultaneously on the disk sur- 
face. Different types of cascades may then contribute to the 
observed 7-ray spectrum from a specific blazar, and to the 
neutrino flux. If multiple reconnection sites do contribute, 
then the neutrino fluxes predicted for a single site should be 
considered as upper limits. However, a flare being produced 



by a single reconnection region seems to be equally probable 
given the likely short duration of each reconnection region 

(~ kac/c). 

It is expected that as a result of magnetic reconnec- 
tion a large amount of energetic plasma may enter the jet. 
Therefore, the cascade processes taking place above the re- 
connection region should also occur in the jet plasma which 
moves with mildly relativistic speeds. Hence, flares of the 
type discussed above may be accompanied by additional ra- 
diation processes caused by this energetic plasma. Recent 
multiwavelength observations of Mrk 501 argue that high 
energy processes in blazars are complicated (see e.g. Pian et 
al. 1998, Protheroe et al. 1998), and more than one different 
mechanism may be responsible for producing the observed 
photon spectrum. 

The motion of plasma, and additionally the longitudi- 
nal component of the magnetic field in the jet (we do not 
include its influence on the final 7-ray spectrum) could pro- 
vide coUimation of the produced radiation in the direction 
perpendicular to the disk surface. 

In conclusion, this model predicts that two types of high 
energy flares may occur if the particles (hadrons, leptons) 
are accelerated to high energies in reconnection regions on 
the accretion disk surface. Pure 7-ray flares can be caused 
either by secondary pairs produced during acceleration of 
protons which collide with the soft disk photons, or by pri- 
mary leptons accelerated in the reconnection region. In this 
latter case, the 7-ray spectrum can extend up to TeV ener- 
gies if the disk temperature is low enough. For large, lumi- 
nous disks, the magnetic field energy can be efficiently trans- 
ferred to relativistic protons which can initiate 7-ray and 
neutrino flares via production of pious in collisions with soft 
disk photons. Therefore, we suggest that blazars of FSRQ 
type, for which 7-ray spectra show some evidence of a cut- 
offs at higher energies (Pohl et al. 1997) and features of 
large accretion disks (emission lines), may become sources 
of GeV 7-rays and ~ lO^'^ eV neutrinos. In contrast, BL Lac 
objects, with spectra extending to TeV energies, should not 
emit high energy neutrinos according to the present model. 
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